Asymptotic aging in structural glasses 
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Using a non-local Monte Carlo algorithm, we study the aging of a fragile glass, being able to 
follow it up to equilibibrium down to 0.89Tmc (Tmc is the Mode-Coupling temperature) and up to 
unprecedentedly large waiting times at lower temperatures. We show that the fluctuation-dissipation 
ratio is independent of the dynamics chosen and is compatible with a phase transition, and that the 
scaling behaviour of the aging part of the correlation supports the full-aging scenario. 
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Aging is found in many complex systems out of equi- 
librium, like supercooled liquids 1^, polymers col- 
loids i3|, or spin-glasses j^, and understanding it is a 
necessary step towards a unified description of such sys- 
tems 1^ 0. After a short transient since preparation, 
a state is reached in which one-time observables (e.g. 
energy, enthalpy) vary extremely slowly, while two-time 
quantities (correlations, susceptibilities) strongly depend 
on the age (or waiting time t^, i.e. the time elapsed 
since preparation) of the system as well as on frequency 
Lj (or the measurement time t). Despite recent efforts, 
our knowledge of aging of real materials is scant in the 
theoretically important regime of large and small fre- 
quency, where universal features should show Two 
issues still open are the scaling of correlations and the be- 
havior of the fluctuation-dissipation ratio. 

Consider observables A and B {B couples to an ex- 
ternal field h). The susceptibility x (i-6- the time in- 
tegral of the linear response R{tw,t + tw) = 5{A{t + 
tw)) / 6h{tuj)\h=f)) and the correlation function C{tw,t + 
tyj) = {A{t + tw)B{tw)) are expected to be of the form 

c(i^,t» + = CstW + agf^^%i^V (1) 

where g(t) is a monotonic function acting as an 'effec- 
tive' correlation time, and Cag describes the aging of the 
system 0. Most published studies focus on the scaling 
properties of Cag: it is generally a function oit/f^, but 
there is a lack of universality in the values of the expo- 
nent /i, embarrassing in view of the claimed equivalence 
of complex systems. Full aging (/x = 1) has been clearly 
observed so far only in spin-glasses 0. For colloids, both 
superagmg and full aging has been reported 

[inj l. Polymers show rather suhaging (/x < 1) [2, 
as has also been observed in simple liquids [l^l- How- 
ever, the values quoted often correspond to different time 
regimes, and the regime where 00 with t/tw fixed 

has not been carefully studied (except for spin glasses). 
For example in glycerol 0| full aging has not been seen 



either close to the glass temperature, Tg (almost at equi- 
librium) or at lower temperatures T. In both regimes the 
explored frequencies were much larger than 

Aging is also characterized by a non-trivial behavior of 
the fluctuation-dissipation ratio (FDR), namely 

X{t^, t + tw) = ^-^(^""^ + ^'"') ^2) 

The fluctuation-dissipation theorem (FDT) states that 
X — \ 'm thermodynamic equilibrium, but this need not 
be so during aging, and FDT violations (i.e. X ^ \) 
are observed. Experiments [l^ Il5| . mean-fleld results 
[lit and simulations [13, E3 suggest that the FDR de- 
pends on time only through the correlation function, i.e. 
X — X\C{tw,tw + 1)]. In structural glasses, in which we 
concentrate from now on, simulations also show that at 
fixed tw, X takes essentially two values: X{C) — 1 for 
C greater than some q^j^{T) (called Edwards- Anderson 
parameter) and X{C) = a;(tu,) < 1 for C < qEA{T). 
Since T/X can be interpreted as an effective tempera- 
ture Teff 01 > it seems that FDT violations in structural 
glasses can be characterized by a single time-dependent 
Tcs{tw) = T/x{tu,), related to the slowest degrees of free- 
dom. This lacks experimental confirmation. (Note that 
other definitions of effective temperautres have been ex- 
plored 0,|l3)- Also open is the issue of the behavior 
of Teff (ill)) as tw 00 (numerical data available cover 
only very short waiting times in the sense that one-time 
quantities are still fastly evolving 0,0|), of great the- 
oretical interest because it is related to the possible ther- 
modynamic meaning of T^g [l^. 

In this paper we study the aging d yna mics down to 
0.53 Tmc [Tmc is the Mode-Coupling |23| temperature, 
below which dynamics slows down dramatically), reach- 
ing very large waiting times. This can be achieved 
through the use of a non-local algorithm (SMC [1^), 
which greatly accelerates the dynamics. We reach an 
asymptotic regime where the correlation function shows 
full aging within errors (supporting the analogy with spin 
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glasses [8j), and where FDT violations are independent 
of the dynamics and of the age of the system. 

We have simulated the soft-sphere binary mixture 
(pair potential Vab{t) = {'^ab/tY'^, diameter ratio 1.2), 
a simple fragile glassformer, using a non-local Metropo- 
lis Monte Carlo algorithm (hereafter SMC) Hil which 
adds swap moves (with probability p) to standard local 
Monte Carlo (LMC). Although swap acceptance is very 
low (« 3 • 10^'^) the equilibration time is considerably 
shortened; e.g. at 0.89 Tmc extrapolations estimate it to 
be three orders of magnitude larger for LMC than for 
SMC (note that other non-local algorithms have proved 
useful in simulations of structural glasses 0|). We used 
the following protocol: Starting from a random config- 
uration, a system of iV = 2048 particles was instanta- 
neously quenched to the final temperature T, and let 
evolve for steps. This preparation was done with the 
SMC algorithm with p = 0.1, which gives the faster equi- 
libration for this system size. After t^, the correlation 
and response functions in the presence of an external field 
h were computed, mostly in SMC runs with p = 0.1, but 
also in LMC and SMC runs with different p in order to 
assess the dependence of the results on the dynamics. 
Due to the swap moves, particle diffusion is not a conve- 
nient observable. Instead, we divided the simulation box 
in Nc cubic subcells and considered the quantity 

1 

= (3) 

a=l 

where = ±1 randomly and Ua is the occupation num- 
ber of subcell a. The side of the subcells was about 
0.35 a AA so that essentially = 0,1. Note that swap 
moves do not change A(t). To measure response, a term 
XNA was added to the Hamiltonian, with A = hkBT (h is 
dimensionless). We considered the correlation C{tw,tw + 
t) = {N A{tw)A{tw + t)), where (. . .) means average over 
both thermal histories and the Cq, together with the in- 
tegrated response kBTxit^.tw -\-t) = {A{tw + i))/h |23 |. 

With SMC we can equilibrate the system down to 
T = 0.89 Tmc- The correlation C{tw,tw +t) shows aging 
up to — 10^, but does not change between = 10^ 
and 10^, which is approximately the region where the 
energy reaches a stationary value (Fig. QJ. We conser- 
vatively estimate the autocorrelation time as the time 
T needed for C to reach the asymptotic value N /Nc 
(~ 0.04), obtaining r = 2 x 10^, much smaller than 
10^ (the total lenght of the simulation). Hence we claim 
that the system has equilibrated, which is further con- 
firmed by the fact that the FDT holds. In contrast, well 
below 0.89 Tmc the system is out of equilibrium up to 
in, = 2x lO'' (our largest observational time). A stretched 
exponential fit of the equilibrium correlation in the late 
Qf-relaxation regime yields a stretching exponent (3 ^ 0.3. 
The equilibrium LMC correlation function does not de- 
cay to N /Nc within the simulated times, hence it is still 
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FIG. 1: Correlation function 0(1^,1^ -\- t) vs. t for T = 
0.89Tmc at t„ = 10^, 10'*, 10^ 10^ (24 samples). Bottom, left: 
integrated response Tx vs. correlation function C at t„ = 10^ . 
Top, right: Energy per particle E vs. t during a SMC quench 
with p — 0.1. Error bars are of the order of point size. 



an open point whether SMC changes the shape of the 
correlations in equilibrium, or whether the two dynamics 
are related by a simple rescaling of time. 

We first address the issue of the scaling of the correla- 
tion during aging at T = 0.53 Tmc (in general far below 
Tg, e.g. for glycerol this corresponds to T 140 K, while 
Tg - 190 K). With SMC we find (Fig. 13) that the corre- 
lations for t^, = 5 X 10^, 5 X 10^ can be made to collapse 
by plotting them as a function of t/t!/^, with /i = 1.05(6), 
compatible with full aging. The collapse applies to the 
aging part (Cag, eq. ^1, which dominates the correlation 
for t/tw > 0.1 Jujtw < 10), as has also been observed 
in spin glasses The two shortest i^'s (inset) can in- 
stead be scaled with fi ~ 0.85. The same value (within 
errors) was found in molecular dynamics simulations of 
the Lennard- Jones binary mixture .12] . so we argue that 
the accelerated dynamics does not affect the scaling. If 
one insists on scaling all curves, it can be done reason- 
ably well using p, 0.9, though this is likely an artifact 
of mixing two different regimes. The relevant point is 
that ~ 1 is seen clearly only for 1^ ^ I and in the 
t ^ tyj region, which is where it is expected to hold 9] , if 
structural glasses share the dynamic properties of spin- 
glasses 0. The failure of full aging for t/tu, ^ 1 is hence 
in agreement with dielectric susceptibility measurements 
in glycerol jl^ . We are not aware of experimental stud- 
ies in the conditions where we find full aging, but such 
measurements are clearly needed. 

A second important result is that although the sus- 
ceptibility and correlation are affected by the choice of 
dynamics, the FDR is not. In fact, Fig.|21shows the ratio 
Toff/T (i.e. the inverse of the FDR) at T = 0.89 Tmc dur- 
ing aging and up to equilibration for both SMC and LMC 
algorithms, obtained measuring the FDR in simulations 
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FIG. 2: C vs. ^ = 1.05(6), for U = 5 x 10^ 5 x lO'^ 

at T = 0.53 Tmc from SMC runs (24 samples). Inset: C vs. 
t/tS;*'^ for = 5 X 10^ 10* (48 samples). 
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FIG. 3: Tcff /r vs. U for SMC and LMC runs at T = 0.89 Tmc 
(16 samples), obtained by a linear fit of the points of the para- 
metric Tx vs. C plots deviating from the FDT line. Errors 
were estimated with the jacknife method 25]. Inset: Tx vs.C 
fortu, = 10* atr = 0.53 Tmc forp = (LMC) andp = 0.1,0.3 
(SMC), N = 20000 (8 samples). 



that used configurations taken along the SMC quench as 
a starting point. After a short transient 10^ steps) 
the FDRs become indistinguishable within errors. At 
T — 0.53 Tmc ^^nd with LMC, we can reach the region 
of FDT violations only for = 10^, so we look at the 
FDR at fixed for LMC and SMC with p = 0.1 and 
0.3, obtaining a good agreement (Fig.|3| inset). 

Finally, we investigate the FDR for large times at 
T = 0.53 Tmc- In Fig. H we plot T^s at = 5 x 10^, 
10^, 5 x 10^ and 5 x 10^ as a function of the instanta- 
neous inherent structure (IS) energy Eis{tw)- We also 
plot Tefi computed according to the IS approach [l8| . 

= dTj/df, where S(/) is the logarithm of the num- 
ber of IS with free-energy /, and d'S/df is obtained as in 



ref. Hi). This idea (which makes no prediction about the 
tuj — > oo limit of Toff) had previously been confirmed only 
in the very early aging regime by molecular dynamic sim- 
ulations 18]. Our results show a reasonable agreement 
even at quite large times. 

The limiting value of Toff as i^; — * oo is of great the- 
oretical interest. If the system eventually equilibrates, 
then Toff T, as we have found for T = 0.89 Tmc- Ap- 
proaches that consider aging a result of critical slowing 
down due to the proximity of a critical point which is 
never reached (beacuse it is located at T = "2^, or 
because of the impossibility to establish a "liquid" long 
range order [i^]) predict this to be the case for all tem- 
peratures. A different view relates the asymptotic value 
of the FDR to a thermodynamic transition described by 
replica symmetry breaking |2^ . Above the transition, 
X{C) is predicted to reach slowly the equilibrium value 
1 (so Toff — > T 0), while below the FDR should remain 
non trivial and Toff tend to a constant > T, since the sys- 
tem never equilibrates. In this scenario the asymptotic 
FDR is claimed to classify complex systems in univer- 
sality classes 0, A third possibility is that FDT 
violations are due to nucleation and slow growth of the 
crystal phase , in which case at long times one expects 
the coarsening regime to be reached, and so Toff — > oo. 

Our results for 0.53 Tmc do not seem to support this 
last possibility. The data are instead compatible with the 
presence of a thermodynamic replica symmetry break- 
ing (RSB) transition [23, since FDR does not seem to 
change between t^, = 5 x 10^ and = 5 x 10^ {Ejs 
are respectively 1.691 and 1.671). Note that this is the 
same regime where the system displays full aging. It 
cannot be excluded that Tes T, but it looks less likely 
if we note that extrapolating i?is(tiu) to i^, — > oo with 
a power-law gives an asymptotic Eis — 1.642. In the 
first approximation the RSB approach predicts that Toff 
equals the transition temperature, which unfortunately 
has been only roughly estimated |2^]. We just observe 
that, at the qualitative level, the fact that the measured 
Teff/T in Fig. ^levels off at a value greater than 1 in the 
late aging regime supports the RSB scenario. 

In summary, we have for the first time studied nu- 
merically the late aging regime of a simple glass-forming 
liquid using local and non-local Monte Carlo (SMC) . We 
find that the scaling of the correlation functions and the 
FDR during aging do not depend on the dynamics. This 
is a strong generalization of the previous finding [3(1 ] that 
equilibrium relaxation in the Lennard-Jones mixture is 
qualitatively identical under different local dynamics (ex- 
cept, as here, for very short times). We have found that 
correlation functions in the late aging regime show within 
errors full-aging scaling, suggesting an equivalence be- 
tween the aging dynamics of structural and spin glasses. 
This should be searched experimentally at frequencies 
comparable or shorter than 1/iui- We also measured the 
FDR while taking one-time quantities closer to assimp- 



4 



4.5 
4 
3.5 

h" 2.5 
2 
1.5 
1 





^\ 5x10^ - 
5x10 ' • 


IS 

- ' simulations < 




[11 

[2] 


3.5 








[3] 
[41 




c 






[5] 




0.5 




□ 




■ ■ 


□ ° T 




[6] 










[7] 





1.64 1.66 1.68 1.7 1.72 1.74 1.76 1.78 



FIG. 4: Teft/T for T = 0.53Tmc vs. the the istantaneous IS 
energy Eis{tw) measured in SMC runs. Eis{tm) is found by 
energy minimization starting from 50 instantaneous configu- 
rations within a small time window (<< tm) around t^. We 
also show the T^ff predicted by the IS approach p^. Inset: 
the Tx vs.C plot for two different 



totic values than in previous studies. FDT violations do 
not imply a thermodynamic transition. However, if a 
transition does exist, there should be a correspondence 
between the asympotic Toff, which is accessible to exper- 
iments, and the order parameter, which is not |3lj |. The 
FDRs measured in experiments [ij, [iSj and simulations 
[lil up to now depend strongly on the age of the sys- 
tem, hence their utility in investigating the existence of 
a transition is still an open point. Here, we have been 
able to reach a regime where Toff has no noticeable time 
dependence. Interestingly enough, it coincides with the 
full-aging regime. At the lowest temperature, the Tcff 
measured over a time window of 3 orders of magnitude 
approaches a finite value, different from the equilibrium 
temperature. This seems little compatible with a criti- 
cal slowing down (Tcff — > T) or the growth of a crystal 
phase (Tcff oo) and favors rather the phase transition 
scenario. Our result suggests that the relevant informa- 
tion to an understanding of aging in structural glasees 
has to be looked for in a regime that so far had not been 
investigated, cither in experiments or in simulations. 
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